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Abstract Classical spin systems with nonadditive long-range interactions are 
studied in the microcanonical ensemble. It is expected that the entropy of such 
a system is identical to that of the corresponding mean-field model, which is 
called "exactness of the mean-field theory". It is found out that this expec- 
tation is not necessarily true if the microcanonical ensemble is not equivalent 
to the canonical ensemble in the mean-field model. Moreover, necessary and 
sufficient conditions for exactness of the mean-field theory are obtained. These 
conditions are investigated for two concrete models, the a-Potts model with 
annealed vacancies and the a-Potts model with invisible states. 

Keywords long-range interaction, mean-field theory, ensemble inequivalence 



1 Introduction 

Long-range interactions appear in several physical systems; gravitational sys- 
tems, non-neutral plasma, vortices in two-dimensional fluid, and so on [6,7, 
12,13]. Theoretically, long-range interacting systems (LRIS) show some pecu- 
liar features both in equilibrium and in nonequilibrium situations. Ensemble 
inequivalence is one of those anomalous properties in LRIS, and has been 
studied for a long time [4,6,7,17]. In short-range interacting systems, several 
statistical ensembles, for instance the canonical ensemble and the microcanon- 
ical ensemble, are equivalent in the sense that these thermodynamic potentials 
are related by Legendre transformation [22]. However, in LRIS, two ensembles 
may exhibit different results, and this is called ensemble inequivalence (see 
also [5] for spin glass systems). For instance, it is known that the specific heat 
can be negative in the microcanonical ensemble, although it is always positive 
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in the canonical ensemble [24]. In short-range interacting systems, the nega- 
tive specific heat implies thermodynamical instability of the state. Since the 
equilibrium state must be stable against any thermal fluctuation, the nega- 
tive specific heat in equilibrium is forbidden. It is nonadditivity that allows 
equilibrium states with negative specific heats in LRIS. Nonadditivity means 
that the energy of two isolated subsystems is not equal to the total energy 
when they are combined [13]. Due to the nonadditivity, the usual argument 
on thermodynamic stability cannot be applied, and the negative specific heat 
and the thermodynamic stability can coexist in LRIS. LRIS also show anoma- 
lous features in rclaxational dynamics such as the existence of the long-lived 
metastable states. In order to understand the properties of those nonequilib- 
rium quasi-stationary states, dynamics of LRIS has been actively studied [1, 
6,7,14,21,25]. These peculiar features in LRIS are ultimately related to non- 
additivity. Properties of systems without additivity have been less understood 
compared to those of short-range interacting systems. Therefore, to under- 
stand many body systems with long-range interactions is an important and a 
challenging problem. 

The so-called mean-field models (infinite range models) are extensively in- 
vestigated because they are tractable models which display those unfamiliar 
characteristics of LRIS. It is expected that general LRIS whose interaction de- 
cays as 1 /r" with < a < d in d-dimensional space exhibit similar properties 
to the mean-field models. Actually, many works support the exactness of the 
mean-field theory in LRIS; that is, LRIS have identical equilibrium properties 
to the mean- field models [3,8,10,19]. However, if the system has some con- 
served quantities, the exactness of the mean-field theory does not necessarily 
hold. In [19,20], it was shown that classical spin systems with long-range in- 
teractions have a temperature region where the mean-field theory is not exact 
in the canonical ensemble with a fixed magnetization. 

In this paper, we study classical spin systems with a long-range interaction. 
We consider systems whose magnetization is not fixed, but whose energy is 
conserved, that is, we treat them in the microcanonical ensemble. Barre dis- 
cussed equilibrium states of the Ising model with long-range interactions in the 
microcanonical ensemble, and argued that equilibrium properties of this sys- 
tem are independent of the power of the interaction a as long as < a < d [2] . 
Campa et al. calculated the partition function of the n-vevctor spin models 
and concluded that the exactness of the mean-field theory holds in these sys- 
tems both in the canonical ensemble and in the microcanonical ensemble [9]. 
We will show that these results do not hold in general. In the microcanoni- 
cal ensemble, the exactness of the mean-field theory is not necessarily valid. 
Necessary and sufficient conditions for the exactness of the mean-field theory 
in the microcanonical ensemble will be presented. We will see that these con- 
ditions are fully described by mean- field quantities and one parameter JJmax- 
The parameter [/max characterizes the interaction potential K{r) and is eas- 
ily calculated. Therefore, we can judge the validity of the exactness of the 
mean-field theory only by analyzing the mean-field model. 
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2 Model 

We consider a system described by the following Hamiltonian: 

N , N 

2 



N N 

H = --J2 Kin - rj)^lJ^j - /I • E (1) 



We consider the d-dimensional cubic lattice, and each lattice point at rj € 
[1, L]'*nZ'* is associated with a spin variable crj. The number of lattice points (or 
spins) is denoted by = L"^. Spin variable "o^ is a general g-componcnt vector 
whose norm is assumed to be finite. We can easily generalize the interaction 
potential K{r) to a matrix form, J2a b ^^abcrfcr^, where a and b are components 
of o^, but we restrict the Hamiltonian to the simple form (1) in order to 
make the discussion simple. The coupling constant matrix J is assumed to be 
symmetric and positive-semidefinite. The external field /i is a g-dimensional 
vector. The interaction potential is long-ranged. We refer to LRIS as systems 
with the following interaction potentials (i) or (ii): 

(i) non-additive limit: the Kac potential is given by 

K{r) = 7'^?!>(7r) > 0, / (l){r)d^r < +oo. 

Jr<fL 

We take the limit 7 — >■ 0, L — >■ oo with 7L fixed {L is the linear dimension of 
the system). Additionally we assume the existence of a decreasing function 
tp{r) such that |V(/)(r)| < ip{r) for all r < jL and there exists a positive 
constant C such that Jg^rK-yL'^^i'''^^''^''^ < C/5 for all 5 > 0. Existence of 
the function ijj restricts the sudden change and the rapid oscillation of 
the potential (j), and the condition ig^j.^^[^'^{i')d'^T < C/5 means that 
ipir) < l/r''+^ in short distance. These conditions are necessary to derive 
Eqs. (9) and (10) discussed later. Note that it includes power-law potentials 
Kir) ~ 1/r", < a < d. 

(ii) van der Waals limit: The potential K{r) is given by 

K{r) = 7'^?!>(7r) > 0, J (l){r)d^r < +00, 

and we take the van der Waals limit, that is 7 — >■ after L — >■ 00. In the 
van der Waals limit, the thermodynamic potential becomes independent 
of the original interaction form (p. Actually, in the van der Waals limit, 
the model is exactly described by the mean-field theory with the Maxwell 
construction [16]. 

An important difference between (i) non-additive limit and (ii) van der 
Waals limit is that the system is nonadditive in (i) but it is additive in (ii). 
The interaction range is roughly given by 7"^. In the non-additive limit, 7I/ 
is fixed at a constant, and therefore the interaction range is comparable with 
the system size. In this case, the system is nonadditive. On the other hand, in 
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the van dcr Waals limit, 7 is taken after the hmit of L — > 00. This means 
that the interaction range is much longer than the microscopic length like the 
lattice spacing but much shorter than the system size. In this case, the system 
is additive and docs not display anomalies such as negative specific heats. 

In the case of (i), the interaction potential depends on the system size L. 
Therefore, we consider the limit of large system size and the limit of weak 
coupling simultaneously in a proper manner. This scaling procedure is called 
Kac's prescription [15]. We normalize the interaction potential as 

AT 

Y^K{r,) = l. (2) 

i 

We can safely take the thermodynamic limit under this normalization because 
the system has extcnsivity (but the system is still no n- additive). 
We refer to the mean- field model corresponding to (1) as 

Hmf = J2 '^l^'^i -lt-J2^i = ^ (—J^'^J^ -t 

= Numf{^), (3) 

where the vector r?t = jfJ2i'^i the magnetization and Umf("^) is the 
energy per spin in the mean-field model. 

It has been revealed that the mean-field thc;ory is exact for the canoni- 
cal ensemble without fixing the magnetization but it is not for the canoni- 
cal ensemble with restriction of a magnetization value in LRIS with periodic 
boundary conditions [19]. Thus it depends on the specific ensemble whether 
the mean-field theory is exact or not in LRIS. In the present paper, we ex- 
amine whether the exactness of the mean-field theory holds or not for the 
microcanonical ensemble with an unrestricted magnetization value in classical 
spin systems. 



3 Microcanonical entropy 

We consider a spin system with a fixed energy and without fixing the mag- 
netization. Such a situation is described by the microcanonical ensemble with 
an unrestricted magnetization value. 

In the microcanonical ensemble, the natural thermodynamic potential is 
the entropy per spin s{e) which is a function of the energy per spin s. The 
entropy is defined by the Boltzmann formula: 

s{e)=lim ^logW{e), (4) 

where W{e) is the number of states with an energy per spin e. We have chosen 
the temperature unit so that the Boltzmann constant is unity. In systems with 
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short-range interactions, the entropy is obtained by the Legendre transforma- 
tion of the canonical free energy 

Scan(£)=inf[^(£ -/(/?))], (5) 

where /3 = 1/T is the inverse temperature. However, in LRIS, the above for- 
mula does not necessarily hold, s(e) ^ Sca.n{s), due to ensemble inequivalence. 
Therefore, we must calculate the entropy directly from the microcanonical 
ensemble. Later, we will distinguish these two entropies, Eqs. (4) and (5). 
The former is referred to as the microcanonical entropy, or more simply, the 
entropy. The latter is called the canonical entropy. 

The free energy in the canonical ensemble is expressed as 

pm=M[ps-s{s)] = s*{p), (6) 

€ 

where the superscript * means the Legendre transformation. The canonical 
entropy (5) is thus expressed as Scaii(e) = s**{s), namely, it is obtained by 
applying the Legendre transformation twice on s(e). It is well known that 
s**(£) is the concave envelope of s(e) (here the "Legendre transformation" 
of a concave function f{x) is implemented by f*{y) = inf^[a;y — f{x)]). The 
concave envelope of s(e) is defined as the smallest concave function greater 
than or equal to s{e) for all e. Therefore, the following relation 

Scan(£) = S**{S) > s{s) (7) 

follows. The canonical entropy gives the upper bound of the microcanonical 

entropy. Since the mean-field theory is exact in the canonical ensemble, that 
is, f{P) = /mf(/3) in LRIS, Scan(e) = Smf(£) holds. Therefore, s(e) < s^-p{e), 
that is, the microcanonical entropy is bounded above by the concave envelope 
of the mean-field entropy. 

On the other hand, the microcanonical entropy of the mean-field model 
smf(£) gives the lower bound of s{s) in general. To show it, let us consider 
the following thermodynamic process. We consider a Hamiltonian H{\) = 
(1 — A)i?MF + Aifi. Here Hi = — ^ J^ij JijCfiCJj is the Hamiltonian of the system 
of interest (Jjj is arbitrary) and Hmf = '^i'^j ^he Hamiltonian of 

the corresponding mean-field model. The paraniciter J is set as J = Jij /N. 
Initially, we set A = and prepare the microcanonical ensemble of the Hamil- 
tonian Hmf with the energy per spin s, which is the initial state of the process. 
Then, we isolate the system from the environment and change the parameter 
X{t) from A(0) = at time t = to A(t) = 1 at time t = t. Since this pro- 
cess is an adiabatic process, the entropy of the final state is greater than or 
equal to that of the initial state regardless of the protocol \{t) of changing A. 
Namely, s{{H)fin/N) > Smf(£)- The expectation value of the energy of the fi- 
nal state is denoted by (iJ)fin- If we change the parameter A suddenly from to 
1 (r is infinitesimal), the expectation value of the energy is (il)fin = 
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The symbol (•)^^ denotes the microcanonical average under the mean- field 
Hamiltonian Hmf and the energy s per spin. It follows that 



ij 

= -lEJij{<ri)f^{<r: 

ij 



MF 
3 Is 



It] 
I'd 

= -^m' = {H^^)f^ = Ne, (8) 

where m = (Yin ai)^^ /N. Thus s{e) > Smf(£) is derived. The same argument 
is applicable to more general Hamiltonians (1). 

The coarse graining of the spin configuration is expected to be performed 
without errors in LRIS because short-scale structure of the local magnetiza- 
tion does not contribute to the free energy per spin. In LRIS, indeed, it is 
known that the free energy in the canonical ensemble is given by the following 
expression [19]: 

/(/3) = min^(/3,{^(aj)}) (9) 

with the free energy functional which is an analytic functional of /3 and 
r^{x). The r^{x) is the coarse-grained magnetization at a point x. For the 
Hamiltonian (1), the free energy functional is given by 



J^{l3,{'^{x)}) = -1- [ d'^x I d'^yU{x-y)l^^{x)Jl^{y)-ll ■ [ d'^x^{x) 
-T [ d'^xuj{l^{x)) 
= U[^{x)]-T [ d'^xuj(T^{x)), (10) 

JCd 

where Cd is the d-dimensional unit cube. The scaled coordinate a; G is 
related to the original position vector by x = Vi/L. The function uj(j^) is 
defined as 

w(rrt) = lim log ( the number of states with a fixed magnetization 

1 ^ \ 



(11) 



which is referred to as the "rate function" in this paper. The scaled potential 
U{x) is given by 

U{x) = lim L'^K{Lx). (12) 

L^oo 
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In the case of (i) non-additive limit, because of 7L — 1, the scaled potential 
is nothing but ^(a;). In the case of (ii) van der Waals limit, U{x) = 6{x). The 
free energy can be expressed by the free energy functional in both cases, but 
the interaction potential U{x) is non-local even after the coarse-graining in 
the non-additive limit. 

In order to derive Eq. (9) in the non-additive limit, the conditions on the 
non-additive limit discussed in Sec. 2, e.g. the presence of ^Pir), is necessary. 
In this paper, this derivation is not discussed. For the derivation, see the 
reference [20]. 

The minimum value of the free energy functional for r^(x) gives the free 
energy in the canonical ensemble. If the magnetization is uniform, r^[x) = r?t, 
the free energy functional is the same as the mean-field free energy with a 
fixed magnetization: 

/MF(/3,7?t) = J-(/3,{77t(a;) = 77l}). 

Namely, if the magnetization is uniform in the space, this state is described 
by the mean-field theory. On the other hand, if inhomogeneity appears in the 
equilibrium state, it is a sign of violation of the exactness of the mean-field 
theory. 

In LRIS, the original Hamiltonian (1) defined on a lattice can be ap- 
proximated by the "energy functional" U[r^{x)]. Indeed, we can show that 
H = NU[r^{x)] + o{N) for LRIS (see [20]). Therefore, the microcanonical 
entropy is expressed as 

s(e) = sup j / u){r^{x))d'^x : e = U\rk{x)] \ . (13) 

The notation sup[^ : B] means that the maximum value of A under the condi- 
tion B. The entropy s{e) is not obtained by the Legendre transformation (5), 
but a similar formula 



s{e) = sup inf {/3 [s - HP, {l^ix)})]} 
holds. By substituting Eq. (10) into Eq. (14), we obtain 



(14) 



sie) = sup inf 

55t(x)/3>0 



'{e-U[fA{x)])+ [ io(^{x))d 



max 



sup inf 



sup inf 



/3(£-W[^(x)])-F / Lo{'^{x))d'^x:e>U[^{x 
{e-U[^{x)])+ [ uji'^ixjjd^x : e <U[^{x)] 



(15) 

In the last equality, we split the configurations {rrl(x)} into those of e > 
U\rk{x)] and e < U\r^{x)] In the former, the minimum for ^ is realized at 



8 



Takashi Mori 



/3 = 0. In the latter, there exists no minimum point for /3. Therefore, it is 
confirmed that the expression (14) is equivalent to 

s(e) = sup j / co{r^{x))d'^x : e > U[r^{x)] \ . (16) 

The expression (16) is equivalent to Eq. (13) in the thermodynamic limit as 
long as negative temperatures arc ruled out. 

If we exchange sup and inf in Eq. (14), it agrees to the canonical entropy 
(see Eq. (5)), 

inf sup {/3 [e - HP. {^{^)})]} = inf [/?(£ - /(/3))] = s^^e)- 

Again the inequality s{e) < Scan(£) is obtained immediately from (14) because 
supinf(-) < infsup(-). On the other hand, if we restrict r^{x) to uniform 
functions of x in Eq. (13), 

s(e) > sup 

which was already derived by the thermodynamic argument. 

In conclusion, we found that the entropy satisfies the inequality 

smf(£) < s(e) < s12f{£)- (18) 

If the mean-field entropy is a concave function of e, it turns out from Eq. (18) 
that the exactness of the mean-field theory holds also for the microcanonical 
ensemble, s{s) = smf(£)- In other words, if the microcanonical ensemble is 
equivalent to the canonical ensemble in the mean-field model, then these two 
ensembles remain equivalent even if the interaction potential is replaced by a 
slowly decaying function. In Ref. [2], the long-range Ising model was studied 
in the microcanonical ensemble and it was shown that the exactness of the 
mean-field theory holds in this model if the ferromagnetic interaction decays 
slower than 1/r'^. Campa et al. [9] also obtained the same result for a family of 
n-vector spin models. It should be noted that the mean-field versions of these 
models satisfy smf(£) = Smf(^)> ^nd the inequality (18) explains the results 
of Ref. [2,9]. 

However, it is known that the microcanonical entropy is not necessarily 
concave in LRIS; negative specific heats are observed in some models. In those 
models, we cannot conclude from Eq. (18) that the exactness of the mean- field 
theory in LRIS holds in the microcanonical ensemble. 
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: e : 



h ■ rrt 



smf(£), 



(17) 



Remark on the van der Waals limit. 

It is verified from Eqs. (13) and (19) that s{e) = Smf(^) in the van der Waals 
limit. Because U{x) = 5{x) in the van der Waals limit, the energy functional 
is of the form, 

{x)] = / UMF{j^{x))df'X. (19) 
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Substituting Eq. (19) to Eq. (13), we obtain 



s{e) = sup < / (jj{rk{x))df'x : s — UMF{'rk{x))d'^x 

T^{x) UCd JCa 



= sup < / sup [u)(rk{x)) : UMFij^ix)) = £{x)] df'x : I £{x)d'^x = e 
/ SMFis{x))d'^x : / s{x)d'^x = s 



= sup 

e{x) UCd 

Last equality is derived by showing ttiat ttie function 



(20) 



sup 

s{x) 



/ SMFis{x))d'''x : / e{x)d'''x = s 

JCd JCd 



is a concave function of e (Proof of the concavity is essentially the same as 
Appendix. C of [20].). The last equality of Eq. (20) is obtained by combining 
this fact with the inequality smf(£) < s{e) < sl^p{e), since Smf(^) the 
smallest concave function greater than or equal to smf(£)- 

Equation (20) means that the microcanonical entropy is obtained by the 
mean- field theory with the Maxwell construction [7]. This result was well 
known one obtained by Lebowitz and Penrose [16]. In later sections, we focus 
on the case of (ii) non-additive limit. 



4 Local stability of the uniform solutions 



We study the local stability of uniform states by the Lagrange multiplier 
method. Hereafter, we impose the periodic boundary condition. First, for con- 
venience, we rewrite the energy functional by the Fourier modes. 



U[rA{x)] = --^ X] UnT^'^Jrk-n - li ■ rrto, 



(21) 



where 



[ d'^xlfi{x)e^'''"-'', 

/ rf''a;[/(£c) cos(27rn ■ cc). 

JCd 



(22) 



Here, Un is called the interaction eigenvalue. The maximum interaction eigen- 
value except for n is denoted by C/max, that is. 



t^max = max?7„. 



(23) 



This quantity plays important roles later. 
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As the total energy is fixed, we search the extremum of the Lagrange 
function 

L({7?t„},/3) = / uj{^{x))d'^x + p{e-U[^{x)]), (24) 

where /3 is a Lagrange multiplier. Let us consider the uniform states, r^{x) = 
rrt. Then the extreme condition is obtained by performing the functional 
derivative 6/6rA{x) and putting rrt(x) = rrl, 

^ w(^) + /3(J^ + 7^) = (25) 



and 

s = ~r^'^Jr^-li -r^. (26) 

If rrt* denotes the magnetization which satisfies Eqs. (25) and (26), the fol- 
lowing relation holds: 

l3{J^*+t) = ~^oj{Tit*)- (27) 
If we regard rrl* as a function of e, we obtain w(rrt*(e)) = smf(£) and 

Therefore, we have 

^=^^-^mf(.). (29) 

This is nothing but the microcanonical temperature of the mean-field model. 
One of the obtained solutions /3 and rrt* corresponds to an equilibrium state 
of the mean-field model. 

Next, we consider local stability of an equilibrium state rrt*(e) of the mean- 
field model. If there is some small fluctuation which increases the entropy 
\^^df'XijjijK(x)'), it is said that this state is locally unstable. The stability 
problem under constraints is solved by analyzing the sign of the minor deter- 
minant of the bordered Hessian matrix (for example, see [18] and Appendix A). 
It is found that the bordered Hessian matrix is block diagonalized into the 
space with different wavenumbers {±n}, that is, 'H = diag('H±n)nez<*,ni>o- 
The block of n = is 



% = \ \ J I . (30) 

\ drrio driiodriio I 

Here d'^L/{drAodrAo) is a q x q matrix whose matrix elements are given by 
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Since the uniform state is stable in the mean-field model, determinant of 
{—l)'^Ho is positive. The block of ±n 7^ is 

ti±n = r.-^ r,-^ = ^^^+pUnJ. (31) 

omnOni-n ornarn 

Because the bordered Hessian matrix is block diagonalized into spaces with 
different wavenumber ±n, the problem is reduced to an eigenvalue problem 
of each block T-L±n- We can judge the stability of the equilibrium state of the 
mean-field model by analyzing only the sign of the maximum eigenvalue of 
{'H±n}- Since the matrix J is assumed to be positive semidefinite, the maxi- 
mum eigenvalue of RHS of Eq. (31) is an increasing function of Un- Therefore, 
the maximum eigenvalue of 

'Hmax = g-^g-^- + PUmaxJ (32) 

determines the stability. The uniform state is locally unstable and the mean- 
field model becomes non-exact if the maximum eigenvalue of the matrix "Hmax 
is positive. Contrarily, if the maximum eigenvalue of 'Umax is negative, the 
uniform state described by the mean-field model is at least locally stable. This 
is a necessary condition for the exactness of the mean-field theory. 

Finally, we mention the relation between the bordered Hessian matrix and 
the microcanonical entropy. Prom Eqs. (27) and (28), we obtain 



drit*drit* 



-/3MF(Jrrt* + li) 



= -puFJ + {J^*+t){J^*-\-lt)'^^^. (33) 
Therefore, each block of the Hessian matrix (31) is expressed as 

n±n = -(1 - Ur,)pMF{e)J + (34) 

where we defined W = {Jr^* -\- 't){Jr^* -\- 't)'^. 

The first term of RHS of Eq. (34) is a ncgativc-scmidefinite matrix. Be- 
cause the matrix W is positivc-semidefinite, the second term of Eq. (34) is also 
negative-semidefinite if d'^SMF{£)/ds^ < 0. Therefore, the maximum eigen- 
value of (31) must be negative in the energy region of d'^supi^) /de"^ < 0. 
Namely, the uniform state described by the mean-field theory is locally sta- 
ble if the specific heat of this state is positive. On the other hand, if the 
specific heat is negative, that is d'^SMF{£)/d£'^ > 0, the maximum eigenvalue 
may be positive. Actually, in the van der Waals limit {Un = 1 Vn e Z''), 
'H±n = Wd'^SMF{£)/de'^ and all the homogeneous states with negative spe- 
cific heat arc unstable. In the non-additive limit (C/max < 1), two energy regions 
appear in general; a region where homogeneous states described by the mean- 
field theory are locally stable, and the other region where they are locally 
unstable. 
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5 Upper bound of the microcanonical entropy 

In the previous section, we derived a necessary condition for the exactness of 
the mean-field theory. In this section, we investigate a sufficient condition by 
evaluating an upper bound of the microcanonical entropy. 

Prom now on, we derive the following upper bound of the microcanonical 
entropy: 

s(e) < sup 

S'>£ 

where /3conv is the inverse temperature where the convexity of the mean-field 
free energy is lost, 

^conv = sup {/3 > : /mf(/3, T^) = /mfW, ^)V?^} , (36) 

where /mf(/^i ^) is the convex envelope of fMpiP, ^) with respect to rrt. This 
relation will be proved in Sec. 5.2. 

Prom this relation, we conclude the following statement, which will be 
uscfid to check whether .s(;) = ,SMF(ff) holds: // the function in the bracket 
of Eq. (35) takes the maximum value at e' = e, we obtain s(e) < smf(£) as 
an upper bound. In this case, combined with s{s) > smf{£), we find that the 
mean-field theory gives the exact microcanonical entropy, s{s) = SMF(e)- 



-(e - s') + SMF(e') 



(35) 



5.1 Graphical meaning of the upper bound 

It helps our understanding to graphically represent the upper bound (35). In 
Fig. 1(a), we depict a typical shape of the mean- field entropy smf(£)- Here, 
we assume that smf(£) =/= •Smf(s) i^ ^^^^ region ei < e < St, and focus on this 
region because s(e) — smf(^) holds trivially outside of this region due to the 
inequality (18). In Fig. 1 (a), e is the inflection point of Smf(£), and thus 

9SMF(e) _ n I \ ^ Q /-\ _ 5 

— = 0mf{s) < /3mf(£) = 
More generally, we deflne 

P = max /3mf{£)- (37) 

Sl<£<Sh 

Prom Eqs. (35) and (18), if 

§^(£ - e') + SMF{e') < SMF{e) (38) 

for all e' > e, the mean-field theory is exact, s{e) = Smf(£)- Inequality (38) is 
rewritten as 

SMF is') < SMF (e) + §^ {£' - s) . (39) 

If we depict the graph of the RHS of Eq. (39) as a function of e', it is expressed 
as a straight line of slope /3conv/^^max which passes through a point (e, SMF(e))- 
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Therefore, Eq. (39) raeans that if this straight hne is above the graph of smf(£') 
for ah e' > e, then s{e) = smf(£)- 

In Fig. 2, we demonstrate this aspect. Let us imagine that we change J/max 
from zero to unity. When J/max is smaU, the slope /3conv/C^max of the straight 
Une is large. In this case, the straight line is above the graph of the entropy for 
all e' > e. This situation corresponds to (A) of Fig. 2. As J/max increases, the 
slope decreases, and at last the straight line touches the graph of the entropy 
and the RHS of (35) becomes maximum at e' = £o 7^ which is described by 
(B) of Fig. 2. 

If the condition 



> P (40) 



is satisfied, the inequality (39) holds for any e such that £/ < £ < e/i. In this 
case, s(e) — SMF(e) holds for all e. Hence, the inequality (40) gives a sufficient 
condition for the exactness of the mean-field theory. In Sec. 6, we will explicitly 
examine the upper bound of the entropy (35) in concrete models. 




5.2 Derivation of Eq. (35) 



Here, we present the derivation of Eq. (35). Let us start with the expression 
(14). Because supinf(-) < inf sup(-), 

sie) = sup sup < inf [/3(e — J-{(i, {r^{x)})\ : / r^{x)dfx = r?t > 

T{l3,{l^{x)}): ( 7^{x)d'^x = T^ 
Jc 



< sup inf < P [ e — inf 



m.{x) 



sup inf [P{e-f(fiM)], 
tA '3>o 



(41) 
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e e' 

Fig. 2 An illustration of the upper bound of the entropy (35). We consider a straight line 
which passes through a point (e, smfC^)) a-nd has the slope /3conv/C^max. (A) the line with a 
small J7max, and (B) the line with a large {/max- In (B), the straight line touches the graph 
of the entropy at e' = eo- 



where /(/?, nt) is the free energy at a temperature T — (3 ^ with a fixed value 
of the magnetization rrt. In the first equaUty, we used 



sup [(•)] = sup < sup 



(•) : 1^ 



The free energy with a fixed magnetization satisfies the following inequal- 
ity [19,20]: 

^) > /mf(/3, ^) - (7„,ax[/MF(/?C/max, ^) " mAfiU^... ^)]. (42) 

Because inf^>o[-] < info<;3<^^„,„/c/max [']' obtain 



sis) < sup inf 

^ 0</?</3„,„/C/„, 



[/3{e-/MF(/3,^?t) 



+ C/max [/MF(/5C/max, ^ /MF(/3t^max, "^)]}] 



sup inf 

^ O<0</Jeo„v/C/„, 



{/?[£- /mf(/3,^)]}. 



(43) 



We consider the case of £ > umf(^) and the case of e < wmf(^) sep- 
arately. First, we consider the case of £ > umf("^)- Because /mf(/3, = 
WMF (7^) - (l//3)a;(?7t), the minimum is realized at /3 = and we obtain 

sup inf [/3(£ - /mf(/3, "^)) : £ > umfCt^)] 



= sup [^(rrt) : £ > wmf(^)] 
= SMF(e)- 



(44) 
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Next, let us consider the case of e < umf{i^)- In this case, the minimum 
is reahzed at /3 = /3conv/f^max- Therefore, 



sup inf 

St 0</3</3eo„v/C/„ 



(£ - /mf(/3, r^)) ■■ e < UMF(rrt)] 



sup 

55t 



= sup sup 



= sup 

e'>e 



max 

/3c, 



-(e — UMF{r^)) + a;(rrl) : e < umf{j^) 
-(e - e') + : umf(?^) = e' 



-(e - e') + smf(£') 



(45) 



When e' = s, 



-{e - e') + SMF(e') = SMF(e)- 



Therefore, 



sup 

e'>e 



-(e - e') + smf(£') 



> smf(£) 



(46) 



Prom Eq. (43), 



s{e) < max< smf(£), sup 



-{e - e') + SMF(e') 



= sup 



-(£ - e') + SMF(e') 



(47) 



In this way, we obtain Eq. (35). 

Finally we comment on the upper bound (35). If we consider the van der 
Waals limit and put JJmax = 1 in Eq. (35). this upper bound exceeds Smp(£) 
which is a "trivial" upper bound of s(e). This fact indicates that the upper 
bound (35) is not the optimal estimation of the entropy for large f/max- How- 
ever, this upper bound is enough to lead us to conclude that the mean-field 
theory is exact for small but nonzero Umax such that /3conv/C^max > 



6 Application to generalized Potts models 

In this section, we apply our result on the necessary and sufficient conditions 
for the exactness of the mean-field theory to two examples, the a-Potts model 
with annealed vacancies and the a-Potts model with the invisible states. 
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6.1 a-Potts model with annealed vacancies 



The Hamiltonian of the a- Potts model with annealed vacancies is given by 

N N 

Eli 



N N 



where Sj S {0, 1, 2, . . . , g} and rij = \ri — rj\. The constant is the normal- 
ization constant, 

1 

'^a — ]~i 

which is independent of a specific site i because of periodic boundary condi- 
tions. This model is reduced to the standard g-state Potts model for £> = — oo. 

Now, we consider the case oiq = 2. If we consider the following correspon- 
dence ^ 



the Hamiltonian (48) is expressed as 

N N 



^=-^E^^^-^i-^E(^Oo, (50) 

ij *J i 

where ("^i)o = 1 — ( o^i)i — ( o^j)2- Correspondence between Eq. (1) andEq. (50) 
is the following: 

^>'=^- -'=(i;)- 

The magnetization is given by 



N ^ Vto2 

Here, we define mg — \ — m\ — mi = Ss^fl/N. The Hamiltonian of the 
mean-field model and the rate function are, respectively, expressed as 



Hmf=N 

and 



(53) 

(rrt) = - ^ Wa log nia- (54) 



2 

a=0 

The bordered Hessian matrix 'Hmax is given by the following 2x2 matrix 
(see Eq. (32)): 



1 

mo 



Umax = I ^^""^"^ ''"mi) "to I ^ ^gg-j 

- ^ /Sf/max - + ;^)^ 



mo 
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Therefore, the maximum eigenvalue of %„ 
by A(£), is given by 



at the energy e, which is denoted 



A(£) = /3C/^ax + 2 




m. 



*2 



(56) 



where r^*{e) is the equiUbrium magnetization of the mean-field model at the 
energy e. 

Next we derive ficom of this model. We define x(^, rrl) as 

a;(/3, rrl) = maximum eigenvalue of the matrix [fi, rrl). (57) 

Prom the definition, ^conv is expressed as 



-^conv — sup 



13 >Q: infa;(/3, rrt) > 

rft 



(58) 



Therefore, it is necessary to compute inf^ a;(/3, rrt). The explicit form of a;(/3, rrt) 
is 



ar(/3,7^) = -l + ^ 




— + — + — -W — -— +32- (59) 



This takes the infimum at mo = 0, mi = m2 = 1/2, and we have 

2 



Therefore, in this model, 



inf x(/3, rrt) 

nt 



-1 + 



(60) 



(61) 



We consider the system (48) on two-dimensional lattice for a = 1. In this 

case, it is known that C/max ~ 0.3. First, we calculate the maximum of u){7Yi.) 
under the condition e = H/N. It determines the equilibrium order parameter 
rrt*(e), the equilibrium entropy SMF(e) = ujiri^*) and the microcanonical tem- 
perature T = l/{dsMF/ds) = 1//3mf(^) in the mean-field model (a = 0). The 
quantity P is obtained from the graph of /3mf(£), and A(£) is also calculated 
from Eq. (56). We calculated these quantities for various parameters (e, D). 

When D > 0.5, the phase transition does not occur. On the other hand, 
when < 0.38, the canonical and the microcanonical ensembles are equivalent 
and the result is trivial. Therefore, we focus on values of D between these two 
values. 

Figure 3 shows quantities for D = 0.421 as functions of e. For this value 
of D, the second order phase transition occurs at £ w —0.27 (see Fig. 3(a)), 
although the first order phase transition occurs in the canonical ensemble (not 
shown) , which is a result of ensemble inequivalence. The negative specific heat 
is clearly observed in Fig. 3(b) just below the transition energy, which is a 
sign of ensemble inequivalence. Figure 3(b) also shows $ ~ 1/0.256 w 3.906, 
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which is below /3conv/C^max ~ 2/0.3 w 6.7. Therefore, for this value of D, the 
mean-field equilibrium states are true equilibrium states for all e and mean- 
field theory is exact. Figure 3(c) shows the maximum eigenvalue A(e) of the 
bordered Hessian matrix. In this case, A(e) is negative for all e, therefore the 
mean-field equilibrium states are always locally stable. 




Fig. 3 (a) Equilibrium magnetizations against energies in the mean-field model, (b) tem- 
peratures against energies in the mean-field model, (c) the graph of A(e). The parameter D 
is set as D = 0.421. 



As far as we calculated, A(e) is negative for all the parameter region of 
(e, D) when we put ?7max = 0.3. Thus the equilibrium states of the mean-field 
Potts model with annealed vacancies is at least locally stable for a = 1 in two 
dimensions. 

On the other hand, we found that the sufficient condition (40) is violated for 
0.49 < D < 0.5. Therefore, it is concluded that except for 0.49 < D < 0.5, two- 
dimensional a-Potts model with annealed vacancies for a = 1 is equivalent to 
the corresponding mean-field model in the microcanonical ensemble although 
the canonical and the microcanonical ensembles are not equivalent. In the 
region of 0.49 ^ D < 0.5, the necessary condition still holds. Hence, the mean- 
field theory may be exact even in this parameter region. 
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6.2 Long-range Potts model with invisible states 

Next, let us consider the long-range Potts model with the invisible states. The 
mean-field version of this model is given by 

J N q 

Hm¥ = XI ^"-^"i m '^<^'>«' 

i,j a=l 

where the spin variable cji is di & {1,2, ... ,q,q-\- 1, ... ,q-\- R\. When i? = 0, 
this model is identical to the standard Curie- Weiss Potts model, which is the 
mean-field version of the g-state Potts model. The R redundant states are 
called "invisible states", which do not interact with other spins (see [23]). 
The order parameter is r^"^ = (mi, m2, . . . , m^) and nia — '^i^cn.a/N . The 
density of invisible states is given by toq = 1 — X^^^iTOq. We rewrite the 
Hamiltonian (62) as 

J 

HMF = ~N-J2ml. (63) 
The rate function a;(rrl) is given by 

u;(rrt) = — ^^TOalogTOa — mologi?. (64) 

a=0 

The Hessian matrix "Hmax is given by 

(^max)ab = - ( —Sab + — ) + /SJUma^Sab- (65) 

\ma Too/ 

The a- Potts model with invisible states is expressed by the following Hamil- 
tonian, 

^ = -^E^'^--.E'^-- (66) 

i,j a=l 

This is similar to Eq. (48). Indeed, two models (48) and (66) are equivalent in 
the canonical ensemble, which is observed by putting D = TlogR in Eq. (48). 
The only difference is that the external field —D J2i ^si,o in Eq. (48) is replaced 
by the entropic term — TOologi? in Eq. (64). However, if we consider in the 
microcanonical ensemble, these two models are not equivalent as we show 
below. 

Here, we consider the system (66) with various values of R on the two- 
dimensional lattice for a — 1 and q = 2. The analysis is almost the same as 
that in the previous subsection. It is noted that A(£) takes the same form as 
Eq. (56). And /3conv is also the same as that of the mean- field Potts model 
with annealed vacancies, /3conv = 2. 

We calculate A(£) and for various {e,R). For sufhciently large R {R > 
260), A(£) becomes positive for some s, and the equilibrium states of the mean- 
field model are locally unstable in this parameter region. This situation is 
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demonstrated in Fig. 4. This is totally different from the case of the a-Potts 
model with annealed vacancies. On the other hand, from the calculation of /3, 
it is concluded that the mean-field theory is exact in this model at least for 
R < 12. There is large difference between the necessary condition R < 260 and 
the sufficient condition R < 12. It implies that the derived sufficient condition, 
^ < /^conv/fAnax IS too strong and not optimal, as is mentioned in Sec. 5.2. 



0.5 

-0.5 ^ 
-1 ■ 
-1.5 ^ 

-2 - 

-0.5 -0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 

e 

Fig. 4 The graph of A(£) for R = 400 and f/max = 0.3. We can see the energy region of 
A(£) > 0. In this region, mean-field equihbrium states are locally unstable. 

We demonstrate typical configurations of the equilibrium state calculated 
by the microcanonical Monte Carlo method [11] in Fig. 5. Clearly, the mag- 
netization is homogeneous for e = —0.1 and R = 400 (Fig. 5(a)), but it is 
inhomogeneous for e = —0.0125 and R = 400 (Fig. 5(b)). In Fig. 5(a) A(£) is 
negative, and in Fig. 5(b) A(e) is positive. It is clearly observed that the ho- 
mogeneous state is unstable in Fig. 5(b). More detailed Monte Carlo analysis 
will be reported elsewhere. 

7 Conclusion and discussion 

We found that the exactness of the mean-field theory in long-range interacting 
systems holds when the canonical and the microcanonical ensembles are equiv- 
alent in the corresponding mean-field model. However, it does not necessarily 
hold in the microcanonical ensemble when these two ensembles are inequiva- 
lent in the mean-field model. We derived the local stability condition for the 
uniform state and clarified the parameter region where the system exhibits 
the non-mean-field behavior. It gives a necessary condition for the exactness 
of the mean-field theory. We also derived an upper bound of the microcanonical 
entropy from which a sufficient condition for the exactness of the mean-field 
theory is obtained. It is remarkable that these conditions are described only 
by the mean-field quantities and one parameter C/max- Therefore, we can judge 
the validity of the exactness of the mean-field theory only by analyzing the 
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Fig. 5 Typical equilibrium configurations for (a) (e = —0.1,R = 400) and (b) {e = 
—0.0125, i? = 400). The system size is L = 100. Gray and black points correspond to 
Si = 1 and 2, respectively (but there are few black points). The empty site corresponds to 
Si =0. 



mean- field model and calculating Umax- Because the mean- field models are 
numerically tractable or exactly solved, this feature of the derived conditions 
is preferable. Wc examined these conditions for the long-range Potts model 
with annealed vacancies and with invisible states. 

In the previous works [19,20], it was shown that long-range interacting 
systems can have non-mean-field equilibrium states in the canonical ensemble 
with a fixed value of the magnetization. In this case, the violation of the exact- 
ness of the mean-field theory is considered to be related to the first order phase 
transition with varying the magnetic field in the canonical ensemble without 
restriction of the magnetization value. In this work, it has been revealed that 
the violation of the exactness of the mean-field theory in the microcanonical 
ensemble should indicate smf SmF' ensemble inequivalence. In- 

equivalence between the canonical ensemble and the microcanonical ensemble 
implies the presence of the first order phase transition with varying tempera- 
tures, because the first order phase transition in the canonical ensemble makes 
the microcanonical entropy non-concave. 

The above discussion based on the results of this work and the previous 
ones [19,20] indicate an underlying general property; if a long-range interacting 

system exhibits a first order phase transition when an ijitnisive parameter x 
is varied, then this system may have non-mean-field equilibrium states in the 
ensemble with a fixed extensive variable conjugate to x. 

In this way, the result of this work shows nontrivial properties of slowly de- 
caying long-range interacting systems and will promote better understanding 
of them. 
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A Bordered Hessian matrix 

In this Appendix, we briefly explain the method of the bordered Hessian matrix. We want 
to find a maximum of a function f{~^) under the condition g{~^) = 0. The vector 1^ is 
assumed to be a vector with n-components. Wc introduce the Lagrange function with a 
Lagrange multipUer A, 

L(A,^) = /(^)-As(^). (67) 
Candidates of maximum are obtained by 



dL _ dL 
aA ~ °' 9^ 



0. 



(68) 



Solutions of Eq. (68) are denoted by A* and lif*. 

The solution (A*, "of*) does not necessarily give a maximum point. In order to judge 
whether the solution gives the maximum, the method of the bordered Hessian is used. The 
bordered Hessian is defined as the following matrix. 



n = 





dg 



_ dg 
dx'dxi 



1 __S2_ a £ 

\ dx* dx*dx^ 



dg \ 



dx ^ dx 



(69) 



The fe-th order minor determinant is defined as det?^^*^', where 



/ 

dg (Pl 
dx'i dx'i dx'^ 




dx'^dx^ / 



(70) 



Whether the point ( A* , af * ) is maximum can be judged from the sign of the minor deter- 
minants of the bordered Hessian matrix. The point (A* , ) is maximum if (— I)*" det 'H*-''-' > 
for all fc = 2, 3, . . . , n. See [18] for more detail. 
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